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Abstract 

The reflectivity of photonic cluster is important property for nanopho- 
tonic applications especially when one needs to hide the cluster or to 
make it extremely visible. Currently, where are no methods clearly 
describing how to create the photonic cluster with the predefined re- 
flectivity. In this paper several analytical methods are proposed to 
create the photonic cluster with minimal or maximal reflectivity. The 
proposed methods are applicable for the clusters made of small parti- 
cles. 

1 Introduction 

Invisible men, time traveling, and teleportation are the very known and fas- 
cinating subjects from ancient tales and modern fiction books. For centuries 
these subjects were something easy imaginable but "not of this world" . While 
time traveling and teleportation are still far from a realization (for living or- 
ganisms at least) the idea of invisibility slowly began to translate into reality 
PQ-[5]. The super reflection can be considered as antipode to the invisibility 
and it has many practical applications already (in Bragg gratings and road 
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signs, for example). In some sense the invisibility and the super reflection are 
closely related but inverted phenomena and that is why they can be studied 
in parallel. 

The discussion about the invisible objects and nonradiating sources was 
started in the scientific literature many years ago (see for example the works 
[B]-[7] and the informative review [8]). Recently it was suggested that meta- 
materials can be used as building blocks for invisible objects and the key 
component of these materials is actually nanostructured photonic crystal [1]- 
0, [9]- [ID]. While it was demonstrated that some devices show relatively low 
reflectivity in some directions and at some wavelengths [3], [5], the transpar- 
ent algorithm for design of invisible photonic clusters is not proposed yet. It 
is worth to mention several methods using supernatural conditions. One of 
them is the transformation optics method where the coordinates are trans- 
formed in such a way that the object becomes invisible [H]-[13]. The catch 
of the method is that transformation of the coordinates is equivalent to the 
redistribution of the permittivity inside the object and as the result of such 
transformation the object becomes invisible. It is not clear however, how 
such method can be applied to the photonic crystals formed by many dis- 
crete particles. Another method based on so-called negative refractive index 
metamaterials requires supernatural conditions (see, for example, the review 
[33] and the refreshing work [T5] ) and it is hardly feasible [2] . Inverse scatter- 
ing methods [H]-[T7] can be used, in principle, for construction of invisible 
objects, however a robust algorithm is not developed yet. 

In this paper I will study visibility of the cluster made of small dielectric 
particles by using the local perturbation method (LPM) approach [TS]- (20] • I 
will present and discuss several methods to construct invisible and extremely 
visible photonic cluster made of independent and interacting scatterers. 

2 The field scattered by the photonic cluster 
in the LPM approximation 

The theoretical framework I use is presented in many works (see for example 
[T5]-[2n] and references wherein) and it will be only briefly presented here for 
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convenience and consistency. Consider the photonic cluster made of particles 
which characteristic sizes are small compared to the incident wavelength A. 
For definiteness it is assumed that the cluster is positioned at the origin of 
the coordinates. The electric field E propagating in the host medium filled 
with iV small particles is described by the following equation [20] 

(A - V <g> V + k 2 ) E(r) + - V E(r„)( £sCi „ - e )/ n (r - r n ) = S(r), (1) 

£ ° n=0 

where 

, 2tt uj — f 1 inside particle , nS 

k = — = — we , j n (r — r n ) = < . (2) 

Ac { 0, outside particle 

Here A and V are the Laplacian and nabla operators respectively, £g> defines 
tensor product, k = |k| is a wave number in the host medium (the |...| 
brackets denote an absolute value), u is the angular frequency, c is the light 
velocity in vacuum, e sc n and Eq are the permittivity of the n-th particle and 
the medium respectively, /„ is the function describing the shape of the n-th 
scatterer positioned at r n , and S is the field source. The characteristic size 
of n-th scatterer is denoted as L n . 

It should be noted that the equation fll]) is an approximate one and it 
is valid when the small scatterers (kL n 1) are considered. The solution 
of the equation (CQ) is presented in the Appendix A and I will use in the 
following discussion the final result presented by the Eq. ( 1A6I) 

where 

R n = \r-r n \ > L n . (4) 

Here R n is the distance between the observer positioned at r and the n-th 
scatterer placed at r n , V n is the volume of the n-th scatterer. 

In many practical cases, the distance between the cluster and the observer 
is much more larger than the size of the cluster and the inequality |r| 3> 
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max(|r n |) is fulfilled. In this case the scattered field (EJ) can be simplified and 
it can be rewritten in the far zone (k |r| ^> 1) in the following form 

/ 2 ikr N—l 

E sc (r) = -i— (/ - 1 ® l) J] E(r n )K( £sc>n - ^K^", (5) 

n=0 

where 

1 = r/r, r = |r| ^> max(|r n |) and kr 3> 1. 

The formula (jSJ) is the main result of this section and it will be used exten- 
sively in the following discussion. The formula shows that the field scattered 
by the cluster is mainly defined by the weighted fields F>(r n )V n (e sc ^ n — Eq) 
inside the small scatterers and by the phase kl ■ r in the scattering direction 
1. Note, that the interactions between the particles in the cluster are present 
in the fields E(r n ). 

3 The visibility of the photonic cluster 

3.1 Criteria of the invisibility and the super reflectiv- 
ity 

The visibility of any photonic cluster depends on the intensity I sc (t) of the 
field scattered by the cluster (and incident on the detector positioned at r) 
and on the sensitivity level £ of the detector. The sensitivity level of the 
detector is defined here as minimal intensity at which a signal is detected. A 
photonic cluster will be invisible (in the frequency span Au and at the cone 
of the directions Al) when J sc (r) is smaller or equal to £. It is convenient to 
write this statement for the invisible cluster in the following form 

J sc (r) = |E sc (r)| 2 < £, (u e Au and 1 £ Al), (6) 

where £ is the sensitivity level of the receiver. The sign < in the expression 
(jSJ) can be used due to the presence of the noise masking the detecting signal 
when I sc (r) = £. 

To define the super reflectivity one should use the criterium based on the 
cluster's reflectivity rather than on the detector sensitivity. The cluster will 
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be qualified as super reflective when the following relation is true 



7 sc (r) ~ 3, (u G Au and 1 e Al), 



(7) 



where 3 is the maximal possible intensity of the field scattered by the clus- 
ter. For the cluster made of small particles the maximal scattered intensity 
can be estimated with the help of the Eq. (J5J) in the approximation of the 
independent particles (when E(r n ) = E in (r n )/(1 + (e SCtn — eo) r y n /3eo)) and 
without depolarization (whenf I — 1 <g> 1 ) E =/E) and it has the form 



3 



E,;, 



k 4 



16ii 2 elr 2 



N-l 



n=0 



(8) 



Here 



1 + (e sc , n - e )^ n /3e ^ 
where r y n is the parameter defined by the shape of the particles and for 
spheres, for example, it is 



7„ = 1 - fc 1 + i~ kL n , ( kL n < !)• 



(10) 



In the most ideal case one can modify the number of the particles N, their vol- 
umes V n , positions r n , and the permittivities e n to make the cluster invisible 
or extremely visible in accordance with the conditions ([6]) or (J7D respectively. 



3.2 The ways to achieve the invisibility or super re- 
flectivity 

The scattering process involves the inhomogeneous medium (photonic cluster 
in our case) and the incident light. That is why the visibility of the cluster can 
be maximized or minimized by the changing the cluster or the incident light. 
The second possibility is discussed in the section 7 where it is shown how 
the given cluster can be hidden or made super reflective by using artificial 
illumination. In the next section I will discuss approaches to achieve the 
invisibility or the super reflectivity of the cluster by modifying its structure. 
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The photonic clusters can be subdivided in two main categories: com- 
posed of independent particles and made of interacting particles. The clusters 
made of the independent particles will be discussed in the section 5 while the 
clusters composed of the interacting particles will be discussed in the section 
6. 

There are several possibilities to satisfy the conditions (EJ) or (J7J) and to 
achieve the invisibility or the super reflectance of the photonic cluster. One of 
the ways is to solve the inequations (JB]) and (J7j) with respect to the positions 
of the particles r n when the fields E(r n ), the number of the particles N, and 
the properties of the particles {e sc ,n and L n ) are known. (Note that the fields 
E(r n ) can be found analytically by solving the system of 3N linear equations 
of type flA5j) .) This way, though clearly formulated, will involve solution 
of the system of nonlinear equations that is not feasible. Alternatively, the 
positions of the particles in the cluster can be modified till I sc will be smaller 
than £ or of the order of 3. This way was investigated in the work where 
encouraging results were obtained. The method closely related to the linear 
sampling method can be also used. It will be described in the section 6 in 
more detail. 

Consider the cluster positioned far from the observer. In this case the 
formula (J5J) is applicable and the invisibility condition can be presented 
in the following form 



1-1(8)1 K 



4:7ie r 



(11) 



and the condition (j7|) for the super reflectivity is 

N-l 



7-1(8)1 K 



IE,; 



^ V n e n 



n=0 



(12) 



Here the vector K is defined as 

N-l 



K = ^2 E(r n )K(£ SC)ri - e )e~ 



ikhr„ 



(13) 



n=0 



The expressions ( TTTT) and ( IT2|) can be considered as the equations with respect 
to the positions r n , V n , and e scn . The solution of these equations should give 
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the properties of the particles (V n and e sc ^ n ) and their coordinates r n in the 
invisible cluster or the super reflector when the observer is positioned at 
the distance r from the cluster in the direction 1. The equations (fTTj) and 
f fT2|) should be complemented by the additional equations for the fields E(r n ) 
and they should contain at least N equations formulated for N observation 
directions 1. It should be emphasized that the equations are nonlinear with 
respect to the unknowns and they are extremely complicated. While the 
equations (iTTj) and (fl~2j) can be solved numerically in principle, a number of 
serious drawbacks and limitations exist reducing the value of the numerical 
solution. The most serious limitations are the complex values of the solutions 
and the excessive number of the solutions. 

Below, the conditions of invisibility and super reflectance will be discussed 
for the photonic cluster made of small particles. 

4 The zero and the maximal scattering by 
the cluster of independent particles 

Consider now the cluster made of N scatterers placed in such a way that 
interactions between the particles are negligible. The interaction between a 
central particle and M particles surrounding it is negligible when the cube 
of the distance between the particles is much larger than the total volume of 
the particles, i.e. when the following condition holds 

R 3 mn » MK/4. (14) 

For more information one can see the formula ( 1B4I) in the Appendix B and 
the required conditions. In this case the fields E(r n ) inside the scatterers can 
be calculated explicitly by using the Eq. fl A5 j) 

Er„ = — — — -p— , 15 

where 7 n are the parameters taking into account the characteristic size L n 
and the shape of the particles. For example, for spheres and for cubes r y n 
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respectively are 



7r. 



In 



I - fc^ ( 1 + i-kL n 



(16) 



1 - A; 2 L^ I 1.52 + i-kL. 
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where L n is characteristic size of the particle (radius of the sphere or half of 
the cube size). Substituting the fields (ITS"]) into the equations ffTTj) and ( 112"]) 
we have 

jv-i , 



n=0 



N-l 



ifcol-r„ 



n=0 



|E in (0)| 



k 2 



N-l 
n=0 



(17) 



18) 



The equations (IT7|) and (ITS!) show that even for the independent particles the 
cluster's visibility is extremely complex phenomenon depending on all parti- 
cles (via the fields E(r n )) and it is extremely sensitive to the parameters of 
each particle (V n , e sc ^ n , j n , and L n ). The equations (fT7|) and (|T8|) are diffi- 
cult to solve and that is why some simplifications are required. Below I will 
use two important simplifications: the scatterers with identical permittivity 



e sc ) and the long wavelength approximation. 



4.1 The cluster made of independent particles with 
identical permittivity 

Consider the case when the particles in the cluster have the same permittivity 
e sc and the incident field is the plane wave, i. e. E in (r) = Ae* kr . In this 
case the expressions ( ITT]) and f lTB]) can be essentially simplified and they can 
be rewritten in the following ultimate form 

N-l 



n=0 







V„ 



(19) 



where k sc = k — k\ is the scattering vector. The formula ([19]) was obtained in 
the assumption that the sensitivity level is £ and that depolarization is not 
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essential in this case. The Eq. ( II 9p shows that the visibility of the cluster 
made of the independent scatterers is solely defined by the phases k sc -r n and 
by the volumes V n playing the role of weight factors. 

The expression (fl~9"l) is actually system of nonlinear equations for the 
positions r n . When the positions r n are known one can produce invisible 
or super reflective cluster. The downside is that even this simplified system 
(fT9"j) is difficult to solve analytically. 

However, for some specific systems the solutions are clearly visible. Con- 
sider, for example, the cluster with central symmetry. In this case the ex- 
pression (1191) transforms into the following one 



(JV-l)/2 
n=0 





Z^n=0 Vr 



(N—l)/2 TA 

(20) 



One set of solutions of the Eq. (TSUI) is clearly visible and for the invisibility 
it is 

k sc -r n = tt(1/2 + n) (21) 
and for the maximal reflectivity it is 

f 2l7l7T , , 

ksc ' rn = \ (2m + l)vr ' (22) 

where m is an integer. 

It is interesting to note that the maximal scattering happens only when 
the conditions (122]) are fulfilled and the scattering is not affected by the 
"weight factors" V n . On the other hand, for the minimal scattering, the 
weight factors V n are extremely important such that invisibility can happen 
in multiple ways (i. e. other solutions are possible where the weight factors 
V n play active role). 



4.2 The cluster of independent particles in the long 
wavelength approximation 

It is interesting to note that despite the independency (defined by the con- 
dition (1141) the particles in the cluster can be placed such that the distance 
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between the adjacent scatterers will be much smaller than the incident wave- 
length such that the following condition will be satisfied 



|k sc - (r n - r m )| < 1, (23) 

where the particles positioned at r n and r m are the adjacent ones. By using 
the condition (1231) the sum (|T7|) and (|T8|) can be replaced by the following 
integral 

i r 

e(r)e iksc ' r dr 



(24) 



where 



^'"■-^"'■ itM^V ' (25) 

Here d is the average period of the cluster. Note that the integral T depends 
on the shape of the cluster which is not known beforehand and the shape 
should be defined from other considerations (practical or guess, for example). 

Since e(r) vanishes outside of the cluster, the integration in ( 124"1) can be 
extended till infinity and T is actually the Fourier transform of the e(r). In 
this case we can present T as 

r = ^|e(-k sc )|, (26) 

where 

e(k sc ) ee — / e(v)e-^dv (27) 

is the Fourier transform of the contrast function e(r). 

The formula ( 12 6 p shows that when the long wavelength approximation 
is valid, the visibility of the cluster made of the independent scatterers is 
defined by the Fourier transform of the contrast function e(r). This property 
allows to construct the permittivity e sc (r) of the cluster by using the following 
formula 

£se(v) = e ( 1 + ^ttttt) • (28) 

It is worth to note that this result resembles the one presented in [8] for 
the invisible scatterer in the first Born approximation. 
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5 The visibility of the cluster of interacting 
particles: the linear sampling method ap- 
proach 

Consider the cluster of arbitrary form made of interacting particles. Surely, 
the cluster can be surrounded by the sphere of the radius p and we can 
subdivide the sphere into small cells positioned at the points r n (|r n | < p). 
Suppose that we know the field scattered by the cluster at the observation 
points r m such that |r m | = r. In this case one can write the following system 
of equations in respect to the fields E(r n ) and the weighted contrasts p n 



U2 ik r ^-J. 

E sc (r m ) = (I - \ m ® l m ) V E(r„,)/i n e- ifcI —, (29) 

n=0 



N-l 

4ire r 

where m = 1...M (M > N) and 

[i n = V n (£ sc>n - E ), l m = r m /r, \r m \ = r. (30) 

The equations f )29|) can be resolved in respect to the multiplication product 
E(r n )/i n . The next step is to find the fields E(i-j) by using the following 
equations 

Efo) = E in (r i ) + E(r i )// i 7 i + (31) 
— [al - bl nj (9 l nj ) ^ E(r n )p n] -— 

a = l + i/kR-l/k 2 R 2 , (32) 
b = -l-3i/kR + 3/k 2 R 2 , (33) 

where 7^ takes into account the shape of the j-th particle and it is described 
by the formulae above. 

When the field E(i\,') is known, the weighed contrast fij is fount from the 
multiplication product E(r n )/i n . The number of the cells (N) will be larger 
than the number of the particles in the cluster, however the contrast of some 
cells should be close to zero such that shape of the cluster (defined by the 
non zero contrast) will be not spherical but similar to the real one. 
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The only open question is the values of the scattered fields E sc (r m ) since 
we need to have definitive values of them to construct the invisible or the 
super reflective cluster. We note that the amplitudes of the fields E sc (r m ) 
are defined in some sense by the Eqs. (11 II) and (Tl2l while the phases are 
not clearly defined. This can lead to ambiguity which can be reduced by 
imposing additional restrictions (from design or engineering point of view, 
for example). 



6 The visibility of the cluster in artificially 
created incident field 

Some applications require to hide a given cluster or make it clearly visible 
(super reflective, for example). One of the ways to do it is to synthesize the 
incident field E in in such a way that the cluster will change its reflectance 
(at the wavelength A in the direction 1). To do this, the fields E(r ra ) inside 
the particles should be found from the following system of equations 

7,2 ikr N—l 

(l- \ m <g> l m ) Y, E(r n )K(e sc ,n - s )e- m ^ = E sc (r m ), (34) 

n=0 

where \ m and r m are the directions and the points respectively in which the 
cluster should be invisible. The system (IMj) consists of linear equations with 
respect to the unknown fields E(r n ) and it can be easily resolved. When the 
fields E(r n ) are known, the scattered field E sc (r n ) can be found by using the 
formula flH]). The incident fields Ej n (r n ) are found from the following formula 
(see Eq. (]A2j) for the reference) 

E m (r n ) = E(r n )-E sc (r„). (35) 

When the cluster is illuminated by the incident field E in synthesized in accor- 
dance with the expression (135]) . the cluster will be invisible or super reflective 
for the observers positioned in the directions 1^ at the wavelength A. 
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7 Conclusions 



The visibility criteria have been discussed for the photonic cluster made of 
small particles. The conditions of the zero and of the maximal scattering have 
been studied for the clusters made of independent and interacting particles. 
It has been shown that in the long wavelength approximation the visibility 
of the cluster made of the independent particles is governed by the Fourier 
transform of the optical contrast. 

The clear algorithm to construct the invisible or the super reflective pho- 
tonic cluster has been proposed for the clusters made of interacting particles. 

The new method to hide a given photonic cluster or to make it extremely 
visible by creating artificial incident field has been presented. 
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8 Appendix A: The LPM formalism 

The equation flTJ 

k 2 

(A - V <g> V + k 2 ) E(r) + - V E(r„)(e sc .„ - e )/ n (r - r n ) = S(r), (Al) 

is easily solvable with respect to the fields E(r n ) when the positions r n of 
the particles are known (it will invoke solution of 3N linear equations for 
each frequency u). The solution of the equation (LA II) can be presented in 
the following form 

E(r) =E m (r) + E sc (r), (A2) 
where the scattered field E sc is 

E sc (r) = - ( J+Z®Z ) VE(r n )(6 sc ,„-s )$„(r) (A3) 

£ ° \ ^ / n=0 
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and 

Here the incident field E in is created by the source S in the host medium 
and it is not important for our consideration (see for example [22] for more 
details). The tensor / is the 3x3 unitary tensor in polarization space and 
r n is the radius vector of the n-th particle. The field E(r n ) is the field inside 
the n-th particle, f n is the Fourier transform of the function f n , and • defines 
scalar product. Note that the integration in Eq. (IA4|) is over infinite three 
dimensional spaces. 

The equations for the fields E(r n ) are found by substituting r = r n into 
Eq. flA2]) 

E(r n ) = E in (r n ) + - J> sc ,„ - e ) / V 2 ' dqE(r n ). (A5) 

It should be emphasized that the formula (1A2|) is rather general one and 
it describes the field in the medium with photonic cluster of arbitrary form 
made of small particles of arbitrary form. 

The formula for the scattered field (IA3I) can be simplified when the dis- 
tance between the observer and the n-th scatterer (Rn) is large, i.e. when 
Rn 3> L n . In this case the integral $ n can be calculated approximately. Note 
also that the integral $ n can be calculated exactly at least for the spherical 
particles. When R n ^> L n the integration in f]A3|) gives 



E sc (r) = ^[1+ ^-yr-) E V(rn)V n (e sc ,n - s )^-, (A6) 

where 

R n = |r-r n | >L n . (A7) 

Here R n is the distance between the observer positioned at r and the n-th 
scatterer placed at r n , V n is the volume of the n-th scatterer. 

k 2 (~ V<g>V\^ , , , ,e ikRn 



E « (r) = 4^( I+ —)h EMV " (B -" - ea) ^r- (A6) 
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where 



Rr, 



|r - r n | > L n . 



(A7) 



Here R n is the distance between the observer positioned at r and the n-th 
scatterer placed at r n , V n is the volume of the n-th scatterer. 



9 Appendix B: The condition of the indepen- 
dence of the particles in the cluster 

Consider the cluster made of identical small scatterers positioned at r n . The 
cluster consists of the central particle (positioned at ro) and the particles 
surrounding it (M particles are placed at the distance R from the central 
particle). The field inside the central particle positioned at ro has the fol- 
lowing form (see Eq. (1A5[) ) 



E(r ) 



2 a ikR 



k 2 e 



M 

V(e sc - £ ) ^ ( al + b\ 0n 



n=l 



where 



AtxeqR 

E in {r )}/(l + {e sc -e )/3e ) 

a = 1 + i/kR- l/k 2 R 2 , 
b= -l-3i/kR + 3/k 2 R 2 , 



Ion E(r n ) + 



(Bl) 



and 



1. 



Of i 



(ro 



/R, R 



(B2) 



(B3) 



Here V and e sc are the volume and the permittivity of the scatterers respec- 
tively. 

The importance of the formula (1B1I) is that it allows us to estimate the 
distance at which the particles are independent, i.e. when the field inside the 
particle located at r will be independent from the particles located at r n . 
This will happen when the term containing E(r„) in Eq. (IBip will be much 
smaller than the incident field Ej n (ro). The analysis of the Eq. fIBlj) shows 
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that for long wavelengths (when kR <C 1) the particles are independent when 
the following condition is satisfied 

R 3 > MV/4, (B4) 

where M is the number of the particles surrounding the central one. The con- 
dition flB4j) is used in the paper to discriminate the clusters with independent 
and interacting particles. We note that this is over estimated condition since 
it was supposed that all the particles interfere in a constructive way. The 
condition flB4j) shows, for example, that for the cluster with M = 6 surround- 
ing spheres (case of simple cubic lattice) the particles will be independent 
when R > 3.3L. 
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